Short time relaxation of a driven elastic string in a random medium 



o 
o 

(N 
< 



S3 
S3 



S3 
O 

o 



> 

OO 

m 
oo 
O 

O 

» 

03 



S3 
O 

o 



X 
S3 



Alejandro B. Kolton, 1 Alberto Rosso, 2 Ezequiel V. Albano, 3 and Thierry Giamarchi 1 

1 DPMC-MaNEP, University of Geneva, 24 Quai Ernest Ansermet, CH-1211 Geneva 4, Switzerland 
2 LPTMS; CNRS and Universite Pans-Sud, UMR 8626, ORSAY CEDEX 91405, France 
3 INIFTA, Facultad C. Exactas UNLP, CONICET, 
Casilla de Correo 16, Sucursal 4, (1900) La Plata, Argentina 

We study numerically the relaxation of a driven elastic string in a two dimensional pinning land- 
scape. The relaxation of the string, initially flat, is governed by a growing length L(t) separating 
the short steady-state equilibrated lengthscales, from the large lengthscales that keep memory of the 
initial condition. We find a macroscopic short time regime where relaxation is universal, both above 
and below the depinning threshold, different from the one expected for standard critical phenom- 
ena. Below the threshold, the zero temperature relaxation towards the first pinned configuration 
provides a novel, experimentally convenient way to access all the critical exponents of the depinning 
transition independently. 
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The study of the dynamics of elastic interfaces in dis- 
ordered media is relevant for diverse experimental situ- 
ations ranging from magnetioii 2 * 3 . or ferroelectric^ do- 
main walls, contact lines of liquid menisci on a rough 
substrate 6 , to crack propagation^. A fundamental prob- 
lem is the response of these systems to external fields 
which pull the elastic interface with a force /. 

A considerable progress has been done in understand- 
ing the steady state dynamics under the applied force. At 
zero temperature, the system is pinned by the disorder 
and the velocity of the interface remains zero up to a crit- 
ical force f c . Above f c the system undergoes a depinning 
transitior&iSiiiii^ii^ and moves with a non-zero average 
velocity. Fisher first viewed the depinning transition as 
a critical phenomenon^: the driving force / plays the 
role of the control parameter, the mean velocity v is the 
order parameter vanishing at f c with an exponent /?, and 
the divergent correlation length £f ~ (/ — f c )~ v can be 
defined from the velocity-velocity correlation function^. 
The analogy of the depinning transition with standard 
critical phenomena has been however recently challenged. 
By studying the low temperature limit of the steady state 
motion it was found that no divergent lengthscale exists 
below threshold^. 

The non-steady dynamics, although experimentally 
relevant, has received less attention. Recently, Schehr 
and Le DoussaU& have investigated this regime for an 
interface initially flat by analyzing two time correlation 
functions, as / — ► /+. Their functional rcnormaliza- 
tion group calculations show that the transient dynamics 
displays universal behavior. Below threshold, numeri- 
cal studies of the zero temperature relaxation towards 
the pinned state have identified a length scale £f di- 
verging with the exponent v at /,jiLi£. In Ref. [l5| it 
was shown that this length does not affect steady 
state properties, but describes transient processes (de- 
terministic avalanches triggered by thermally activated 
events) during the steady state low temperature motion 
for / < f c , and it is ultimately related to the vanishing 
of the density of metastable states as one approaches f c . 



The / = relaxation towards equilibrium at finite tem- 
perature has been studied both numerically 19,20 - 21 ' 22,23 
and analytically 24 i 25 . 

In the present paper we analyze the transient dynam- 
ics of an interface pulled with a finite force close to f c , 
and show that it is a powerful method to extract the crit- 
ical properties of the depinning transition. Indeed, the 
analysis of the relaxation dynamics has been used exten- 
sively to study equilibrium critical phenomena 26 . The 
basic idea behind this dynamic approach is the existence 
of a growing length L(t). Let us use for simplicity the 
example of the Ising model: the system is prepared in 
the ground state, characterized by a global magnetiza- 
tion m = 1, and at t — it is annealed at a temperature 
T, close to the critical point T c . The global magnetization 
m relaxes to its equilibrium value following a time evolu- 
tion controlled by L(t): for lengths below L(t) the system 
is equilibrated, while for lengths larger that L(t) the sys- 
tem keeps memory of the t = initial condition. After a 
microscopic time, the relaxation is governed by the dy- 
namical exponent z and L(t) ~ t x l z before reaching the 
equilibrium correlation length. In this macroscopic time 
regime scaling arguments lead to a universal behavior for 
the relaxation of the order parameter, although the sys- 
tem is far from its equilibrium stated. Analogously, if 
we assume the presence of such a growing length L(t) in 
the transient regime of an initially flat driven elastic in- 
terface, the scaling form for the relaxation of the velocity 
is given by 



v{t,f)=Cf f,,v F{L(t)ltf) 



(1) 



where the function F(s) ~ s~^l v for small s. When 
s > 1, F(s) ~ const, for / > f c , in order to get the 
steady state velocity v ~ (/ — fc) 13 ] for / < f c , where the 
order parameter is zero at T = 0, F(s) must be modified 
to take into account the exponential decay of the velocity. 
In standard phase transitions, a scaling form equivalent 
to describes, in general, the evolution of the order 
parameter. In this paper we show that the scaling form 
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FIG. 1: (a) Structure factor Sf(q) of the pinned configu- 
ration obtained by relaxing an initially flat line for different 
forces / = 0.325, 2.0, increasing from the bottom to the top 
curve. Inset: collapse using the crossover length £/, and the 
depinning roughness exponent £ = 1.25. (b) Finite size study 
of vs /. The solid line is a fit to ~ (/ c — f)~ v , where 
v — 1.33. The dashed-dotted line indicates the extrapolated 
value of f c . 



(|T|l describes the relaxation near depinning, but in a non- 
standard way. While in standard critical phenomena £/ 
represents the correlation length on each side of the tran- 
sition, for depinning £r represents the steady state cor- 
relation length only above threshold but a purely tran- 
sient correlation length below threshold, absent in the 
steady state geometry of the linei^. With this identifica- 
tion of the relevant lengths the analysis of the transient 
gives access to all the critical exponents of the transi- 
tion. To make the study more transparent we consider 
only the case of strictly zero temperature, in which the 
interface relaxes deterministically towards a final pinned 
configuration without reaching any steady state below 
threshold, thus probing transient deterministic dynamics 
only. We also restrict here to the simple case of a string 
with short range elasticity moving in a two dimensional 
random landscape, but our method naturally applies to 
higher dimensional systems as well. 

The string is described by a single valued function 
u(x,t), which measures its transverse displacement u 
from the x axis at given time t. The equation of mo- 



tion is given by, 

d t u{x,t) 



d 2 x u{x,t) + F p {u,x) + f 



(2) 



where F p (u,x) is the pinning force with correlations 
F p (u, x)F p (u' , x') = A(it — u')S(x — x'), the overline rep- 
resents the average on the disorder realization and A(x) 
is a short-ranged function. In our simulations we have 
set this range equal to 1, analyzed system sizes up to 
L = 2048, and results were averaged from 2000 to 10000 
disorder realizations. For a fixed force / < f c we let 
evolve an initially flat line u(x, 0) = 0, up to its final 
pinned configuration, which can be detected by a very 
efficient algorithm in polynomial time^i. The geometri- 
cal properties of the line can be described using the aver- 
aged structure factor, which for a general configuration 
is defined as, 



Sf («,*) = 



— I dz u(x, t)e lqx 
L Jo 



(3) 



where L is the length of the line, q = 2wn/L, with 
n = 1, . . . , L — 1. For a self affine line we have S(q) ~ 
g _ 0+ 2 C) thus yielding the roughness exponent (. In 
Fig. Ola) we show the structure factor of the pinned line 
Sf(q) = Sf(q, oo ) for different values of / < f c . Two 
regimes can be identified: at small lengthscales, the ge- 
ometry of the line becomes self-affine and it is charac- 
terized by the depinning roughness exponent ^ ss 1.2E 34 . 
At large lengthscales Sf(q) reaches a plateau which rep- 
resents the memory of the initial flat condition. As shown 
by the perfect collapse in the inset, Sf(q) is governed by 
a single length £/, given by the crossover between these 
two regimes. In Fig.^b) we see that the crossover length 
increases with the force, diverging with the depinning ex- 
ponent v at a finite value, identified with the threshold 
f c . Small deviations from this behavior are observed at 
very small forces and whenever the crossover length ap- 
proaches the system size L. Using the finite-size anal- 
ysis of Fig. ^b) we can easily extrapolate the value of 
the critical force for the infinite system and identify the 
crossover length with This length corresponds to the 
size of the minimal string rearrangement needed to reach 
a metastable configuration from the flat one, and its di- 
vergence is due to the vanishing density of metastable 
states approaching f c from below. can be associated 
with the divergent length found in Refs. and with 

the size of the deterministic avalanches in the steady state 
motion 1 -. In the following we show that although does 
not affect steady state properties below thresholcU^ it af- 
fects the transient relaxation. 

Let us now analyze the time evolution of the line to- 
wards the final pinned configuration for / < f c and to 
the sliding steady state for / > f c . For this purpose we 
solve numerically equation by using a second-order 
Runge-Kutta method. In Fig. [21 (a) we show the typical 
evolution of Sf(q, t) for a force / = 1.80. Once again, two 
roughness regimes are observed, one corresponding to the 
memory of the flat initial configuration and the other to 
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FIG. 2: (a) Typical evolution of the structure factor Sf(q,t) 
during the relaxation at / = 1.80, calculated at different 
times, increasing from the bottom to the top curve, obtained 
by using lines of size L = 1024. Inset: collapse using the 
crossover length Lit), (b) Growth of L(t) for / = f c (a), 
/ = 1.75 < f c (□) and / = 2.05 > f c {<>). The lines indicate 
the power law behaviors in the different regimes: L(t) ~ t 2 ^ 3 
(solid), L(t) ~ t 1/2 (dashed), L(t) ~ cte. (dashed-dotted). 
The data for each force has been vertically shifted for clarity. 
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FIG. 3: (a) Relaxation of the velocity v(t, f) for different 
forces around the critical force, from / = 0.5 to / = 2.50. 
The dashed line indicates the power-law behaviour f _/3//l/z , 
where (3, v, z are critical exponents, (b) Collapse of all forces 
using the scaling of Eq.Q based on the existence of a growing 
length L(t) and a second static length associated with the 
density of metastable states. The dashed line indicates the 
power-law behaviour Lit)~^l" . The inset shows the scaling in 
the microscopic time regime, as predicted by the perturbation 
theory (Eq@). 



the depinning roughness = 1.25. As we can see in the 
inset, a growing lengthscale, identified with Lit), governs 
the time evolution. In Fig. [21(b) we show the time evolu- 
tion of L{t) for forces above and below the threshold. We 
can distinguish three regimes for the evolution of L(t). 
After a first microscopic time regime where the line is 
practically flat, we find a macroscopic short time regime 
where the growth of Lit) is controlled by the depinning 
dynamical exponent z ~ 1.5, as Lit) ~ t x / z . This result 
shows that the depinning transition is characterized by 
a universal short time relaxation. The crossover to the 
third regime occurs when Lit) ~ £/, after which we can 
distinguish the relaxation above or below threshold. For 
/ < f c , Lit) saturates to £/, the characteristic length 
of the final pinned configuration. For / > f c , Lit) con- 
tinues to grow as L(t) ~ t 1 / 2 . The thermal dynamical 
exponent z — 2 is produced by the finite velocity, which 
makes disorder act as an effective thermal noise above the 
lengthscale Finite-size effects are expected if one of 



the two important length scales, £/ and L(t), become of 
the order of the system size L. 

In the steady state the observables are controlled by 
the correlation length In particular, the mean veloc- 
ity, order parameter of the depinning transition, is given 
by v ~ £/ '• However, in the transient regime, L(t) 
plays a crucial role in the relaxation of all the observ- 
ables. In particular, in FigEl ( a ) we see that the three 
growth regimes previously found for Lit) correspond to 
three regimes of v(t,f). At the shortest times the line 
moves freely, v ~ /, since pinning forces average to zero 
for a flat line. The first non zero corrections to this free- 
flow behavior can be therefore obtained by using a second 
order perturbation theory in the disorder. This gives, 



A(0) - A(/f) 
/ 



(4) 



The inset of Fig. E^b) shows that (@J holds up to times 
t < rf/f, where rr = 1 is the correlation length of the 
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disorder. The first stage of microscopic relaxation of 
v(t, /) gives thus access to the shape of the disorder corre- 
lator A(it) for u < Tf. The data obtained at longer times 
verify the scaling relation proposed in , as we show in 
Fig. GJb) ■ This collapse shows that the velocity is con- 
trolled by the ratio L(t) /£/. When L(t) <C £/ the growing 
length L(t) is the only relevant scale in the problem and 

replaces in the steady state relation v ~ '. This 

leads to v(t, /) ~ L(t)~^^ u : as it is observed in the initial 
slope of the collapsed curves in Fig. |3Jb), with (3 ~ 0.33. 
Subsequently, for / > f c the scaling function tends to 
a constant while for / < f c has a fast decay. A similar 
relaxation of the velocity was found in other disordered 
models^. 

In conclusion, we have identified the relevant static 
and dynamical lengths, L(t) and controlling the dy- 
namics of relaxation of an elastic line under an applied 
force. We have shown that describes well the uni- 
versal relaxation of the velocity observed for t > 1/f 
and / close to the threshold. In contrast to standard 
critical phenomena, below threshold the static length £t 
does not represent any steady state correlation length. 
The proposed scaling relationship for the time relaxation 
can be used to get all the critical exponents of the de- 
pinning transition independently: we can get u from the 
characteristic length of the pinned configuration below 
threshold (or from the crossover to the thermal growth 
of L(t) above threshold), z from the evolution of Lit), 
and (3 from the relaxation of the velocity. This is best 
accomplished close to the depinning threshold since £/ 
is large, and the macroscopic short time regime longer. 
As for equilibrium critical phenomena, this is a conve- 
nient method for both, numerics and experiments, since 
it avoids the problem of equilibrating a sample in the 
steady state near the depinning threshold. When L{t) 



becomes of the order of the correlation length £f all the 
observables approach their steady state value for / > f c 
or their value at the final pinned configuration for / < f c 
exponentially fast. This study is therefore relevant for ex- 
perimental situations whenever T is low enough to assure 
well separated time scales for deterministic and thermally 
activated motion. 

Finally, we point out that the experimental charac- 
terization of interfaces pinned below f c , prepared un- 
der a unique protocol, should be a convenient tool to 
investigate depinning, since the geometry of these inter- 
faces yields an independent measure of the exponents £ 
and v, as shown in Fig. ^ This feature is important, 
as these exponents are usually related. The so called 
statistical tilt symmetry (STS) is present whenever the 
elastic force is a linear function of the deformation field 
u{x). In general, the STS relation depends on the range 
of the elastic interactions: for a short-ranged elasticity, 
as in , the relation reads v = 1/(2 — Q, while for the 
long-ranged interactions, expected for crack propagation 
in solids or for the contact line of liquids, the relation 
becomes v — 1/(1 — C)- Recent experiments in these 
two systemsL^S gave roughness exponents ( that are sys- 
tematically bigger than the one computed by numerical 
simulations^,. In order to explain these deviations the 
presence of non-linear elastic corrections able to change 
the universal behavior of the depinning transitionSLSLS 
has been invoqued. The measurement of an eventual vi- 
olation of the statistical tilt symmetry relation from the 
study of interfaces pinned below f c could be thus used as 
a novel smoking gun test for this hypothesis. 
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